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The  Linear-combination-of-Gaussian-type-orbitals  (LCGTO)1  local-den¬ 
sity-functional  (LDF)2  ground  state  geometry  of  the  13-atom  iron  cluster  has 
approximate  icosahedral  (I^)  symmetry2  as  does  the  13-atom  nickel  cluster*. 
Zero  and  low  magnetic  moment,  nontransition-metal  elemental  13-atom  clus¬ 
ters  are  also  icosahedral  if  directional  bonding  is  not  too  important2-®.  The 
icosahedral  13-atom  cluster  has  two  symmetry-inequivalent  sets  of  atoms,  a 
unique  central  atom  surrounded  by  twelve  symmetry  equivalent  atoms  each 
occupying  the  corners  of  an  icosahedron.  Thus  either  of  these  two  sets  of 
atoms  can  be  replaced  by  sets  of  atoms  of  a  different  element  without  chang¬ 
ing  the  symmetry  of  the  clusters.  This  work  considers  the  effect  of  the  identity 
of  the  central  atom  on  the  stability  and  magnetic  moment  of  the  cluster  as  a 
whole,  and  how  this  approach  can  be  extended  to  larger  icosahedral  clusters 
that  have  more  symmetry-inequivalent  shells,  which  would  afford  greater  pos¬ 
sibilities  of  alloying  within  the  constraint  of  icosahedral  symmetry.  This  high 
symmetry  can  be  used  to  enhance  the  magnetic  moment  of  these  clusters. 

The  Gaussian  orbital  basis  set  for  iron  was  a  15s/llp/6d  primitive  basis7 
contracted  according  to  an  atomic  calculation  3, 5/2, 3/1, 3,  where  for  each 
orbital  angular  momentum  the  first  number  is  the  number  of  lowest-energy 
atomic  orbitals  (Is  to  3s,  etc.)  and  the  second  number  is  the  number  of  most 
diffuse  Gaussians  used  in  the  cluster  calculations  (5  s-type,  etc.).  Compatible 
primitive  Gaussian  basis  sets  from  a  handbook8  were  augmented  with  a  single 
d  polarization  function9  (/  polarization  functions  on  transition  metal  atoms  in 
compact  structures  lead  to  numerically  singular  overlap  matrices4)  and  every 
exponent  greater  than  10  bohr-2  contracted.  For  sulphur  this  approach  lead 
to  a  12s/9p/ld  basis  orbital  basis  set  which  was  contracted  2, 3/1, 4/0,1.  The 
Perdew-Zunger10  LDF  was  used.  During  the  self-consistent-field  (SCF)  pro¬ 
cess  each  one-electron  energy  level  was  exponentially  broadened  by  0.025  eV 
and  then  occupied  from  the  lowest  energy  up  to  a  single  self-consistent  Fermi 
energy,  c/?,  for  both  spins4. 
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The  number  of  impaired  spins  as  a  function  of  the  distance  between  the 
two  symmetry-inequivalent  atoms  in  the  I*  Fei3  cluster  is  plotted  in  Fig.  1 
and  compaired  with  similar  curves  for  two  less  stable  Fei3  structures.  The 
0h  structure  is  the  cuboctahedron,  which  is  a  compact  portion  of  the  fee 
crystal  lattice.  The  D3A  structure  can  be  obtained  from  the  fee  cluster  by 
rotating  any  triad  of  nearest-neighbor  surface  atoms  by  60°about  their  center. 
This  third  cluster  is  a  compact  portion  of  the  hep  lattice.  Each  surface  atom 
of  the  hep  and  fee  clusters  sees  an  identical  nearest-neighbor  environment. 
Therefore  these  two  less  stable  clusters  would  be  expected  to  have  similar 
properties.  Furthermore,  these  three  cluster  have  the  same  equilibrium  radial 
bond  distance  of  4.4  bohr2.  Nevertheless,  Fig.  1  shows  that  these  three  clusters 
have  quite  different  spin  properties.  For  most  radial  bond  distances  the  cluster 
magnetic  moments  are  related, 


(!) 

This  difference  has  been  attributed  to  different  symmetry-required  degeneracy 
in  these  clusters3.  The  average  degeneracy  of  an  irreduceable  representation 
(IR)  of  the  icosahedral  (I&),  octahedral  (0 a),  and  D3 a  point  groups  are  3.2, 
2,  and  l|,  respectively.  Furthermore,  the  maximal  degeneracy  of  IRs  of  these 
groups  are  5,  3  and  2,  respectively.  The  optimal  spin  for  each  cluster  changes 
quite  abruptly  as  a  function  of  radial  bond  distance  for  these  three  clusters 
because  the  room-temperature  level-broadening  parameter  is  small  compared 
to  the  average  spacing  between  the  one-electron  energy  levels  of  at  least  one 
spin.  For  example  a  jump  of  10  units  between  4.5  and  4.6  bohr  is  required 
for  the  Ift  cluster  because  the  highest  occupied  molecular  orbital  (HOMO) 
and  lowest  unoccupied  orbital  (LUMO)  have  opposite  spin,  h^  and  h„  orbital 
symmetry,  and  reverse  their  order  in  energy'  as  the  radial  bond  distance  changes 
from  4.5  to  4.6  bohr.  Beyond  the  plateau  reached  at  largest  bond  distances 
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in  this  figure,  the  spin-polarized  core  d  electrons  reduce  their  overlap  and 
ferromagnetic  ordering  may  cease. 

Fei3  deviates  slightly  from  perfect  geometry  through  the  Jahn-Teller 
mechanism.  The  HOMO  of  minority  spin  is  a  five-fold  degenerate  hu  orbital 
that  is  occupied  by  only  a  single  electron.  If  the  symmetry  unique  central 
atom  is  removed,  an  icosahedral  shell  of  12  iron  symmetry-equi valent  atoms 
remains.  This  I ^  shell  has  38  unpaired  spins  and  a  radial  bond  distance  of 
4.2  bohr,  yielding  a  radial  contraction  of  0.2  bohr  upon  removal  of  the  central 
atom.  Thus  the  net  spin  of  the  central  atom  is  opposite  that  of  the  shell. 
This  and  other  comparisons  are  made  in  Table  1,  where  the  binding  energy 
is  defined  as  the  total  energy  of  the  13-atom  cluster  minus  the  energy  of  the 
shell  and  an  isolated  central  atom.  These  energies  are  computed  in  (spatial) 
symmetry-restriced  fashion11. 

The  size  of  the  hole  left  inside  a  I/i -symmetric  Fei2  shell  can  be  computed 
two  ways — using  Slater  atomic  radii  or  the  shell  nearest-neighbor  distance. 
The  Slater  atomic  radius12  for  iron  is  1.40A(2.65  bohr),  leaving  a  hole  of 
radius  0.82A(1.55  bohr)  vacant  in  the  center  of  the  Fei2  shell.  The  nearest- 
neighbor  distance  of  two  atoms  in  the  spherical  shell  is  105.15%  of  its  radius. 
Subtracting  half  of  this  nearest-neighbor  distance  from  the  shell  radius  yields 
a  hole  radius  of  l.Q5A(1.99  bohr).  The  small  size  of  this  central  hole  and  the 
fact  that  the  shell  lacks  two  electrons  from  having  completely  filled  or  half- 
filled  shells  suggest  investigating  divalent  atoms  that  are  smaller  than  iron  as 
interesting  candidates  for  filling  that  hole.  Table  1  is  heavily  weighted  with 
such  atoms.  From  this  table,  the  energetically  most  favorable  such  atom  is  zinc, 
the  cluster  containing  which  has  40  unpaired  spins.  Indeed  it  is  the  most  stable 
cluster  found  in  this  study.  The  highest  spin  system  found  is  Fe^S,  but  it  is 
energetically  less  stable  than  Fe^Zn,  although  more  stable  than  the  pure  Feij 
cluster.  If  Fe^O  has  a  slightly  larger  radial  bond  distance  it  too  would  have  44 
unpaired  spins  and  be  stable  against  the  Jahn-Teller  mechanism.  Interesting 
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semiclosed-shell  nondivalent  systems  are  Fe^Al  and  Fe^Mn.  The  former  has 
high  spin,  41  unpaired  spins,  and  the  latter  has  low  spin,  35  unpaired  spins. 


Clearly  these  13-atom  I*  systems  are  rich  in  their  magnetic  properties  and 
thus  could  potentially  be  tailored  for  engineering  applications.  Larger  highly 
symmetric  systems  would  afford  greater  flexibility.  Studies  of  larger  systems, 
specifically  the  Fess  1^  and  Fesi  0&  (bcc)  clusters,  are  under  way. 

With  larger  magnetic  systems  the  problem  of  double  SCF  convergence 
(simultaneously  optimizing  the  one-electron  orbitals  and  the  number  of  un¬ 
paired  spins)  looms  larger.  A  possible  simplification  in  this  effort  is  to  change 
from  using  a  spin  unrestricted  (ULDF)  formulation  of  the  problem  to  a  re¬ 
stricted  open-shell  (ROLDF)  formulation.  Rather  than  direct  optimization 
of  the  ROLDF  problem,  which  in  general  leads  to  off-diagonal  entries  in  the 
one-electron  eigenvalue  matrix,  e,  an  average  of  the  local  potentials  for  both 
spins  weighted  by  the  number  electron  for  each  spin  was  used, 

q&(r)  =  Hv2  +  (2) 

where  the  arrow  subscripts  denote  the  two  spins  and  N  and  V(r)  are  the 
number  of  electrons  and  local  potential.  For  Fei3  the  two  methods  yield  the 
same  ordering  of  the  one-electron  levels  and  thus  the  same  number  of  unpaired 
spins.  For  Ij,  Fess,  with  the  nearest-neighbor  distance  between  shells  fixed 
at  4.4  bohr  (the  optimal  distance  for  Fei3,  but  likely  too  short  for  larger 
systems,  thus  underestimating  their  magnetic  moments)  Figs.  2  and  3  give 
the  ROLDF  and  ULDF  eigenvalue  spectra.  Both  calculations  yield  the  same 
number  of  unpaired  electrons,  64,  but  the  spectra  and  level  ordering  are  slightly 
different.  The  ROLDF  approximation  does  not  significant  improve  the  SCF 
process,  despite  stabilizing  the  orbitals  of  minority  spin  near  the  Fermi  energy. 
Nevertheless,  because  more  than  a  hundred  iterations  are  needed  to  reach 


an  acceptable  level  of  convergence,  the  ROLDF  saves  almost  a  third  of  the 
computer  time. 

Symmetry  can  have  a  significant  effect  on  the  magnetic  moments  of  small 
clusters.  Currently  this  effect  is  most  easily  addressed  theoretically;  thus  it 
is  important  to  extend  theoretical  methods  to  reliably  and  economically  treat 
larger  clusters. 
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TABLE 


Table  1.  Stability,  size,  and  magnetic  moments  for  selected  central  atoms  in  an 
icosahedral  Fe^X  cluster.  The  binding  energy  is  in  hartrees  and  radial  bond 
distance  in  bohr.  An  up  arrow  subscript  on  the  highest  open-shell  spin-orbital 
indicates  that  it  is  occupied  by  majority-spin  electrons. 


Central 

atom 

Binding 

energy 

Radial  bond 

distance 

Unpaired 

electrons 

Open-shell 

electrons 

Open-she 

orbital 

None 

0.000 

4.2 

38 

3 

hg\ 

Be 

0.178 

4.3 

40 

0 

- 

0 

0.118 

4.3 

42 

4 

hg\ 

Mg 

0.010 

4.5 

40 

0 

- 

A1 

0.145 

4.4 

41 

0 

- 

Si 

0.258 

4.4 

40 

0 

- 

S 

0.053 

4.5 

44 

0 

- 

Ti 

0.262 

4.5 

38 

2 

hgl 

Mn 

0.206 

4.5 

35 

0 

- 

Fe 

0.006 

4.4 

34 

1 

Ni 

0.211 

4.4 

38 

3 

kg'] 

Zn 

0.372 

4.4 

40 

0 

- 

Se 

-0.046 

4.6 

38 

3 

h-vi 

Pd 

0.109 

4.6 

38 

3 

h9l 
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FIGURE  CAPTIONS 


Figure  1  The  optimal  number  of  unpaired  spins  as  a  function  of  radial 
bond  distance  for  the  hep  (long-dashed  curve),  fee  (solid  curve), 
and  I/,  (short-dashed  curve)  Fei3  clusters. 

Figure  2  Comparison  between  the  a)  restricted  open-shell  and  b)  unre¬ 
stricted  one-electron  eigenvalue  spectra  for  I&  Fess.  In  each  panel 
the  majority-spin  spectra  are  to  the  right  of  the  minority-spin 
spectra.  Both  methods  give  the  same  number  of  unpaired  elec¬ 
trons,  64,  which  is  fixed  by  the  large  gap  in  the  minority-spin 
spectra  at  the  Fermi  energies. 
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